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Abstract 

The traditional physical description of neutrino flavor conversion in the Sun focuses on 
the notion of resonance. However, the resonance picture is valid only in the limit of small 
mixing angles 9. For large values of 6, the resonance picture leads to seemingly paradoxical 
results. This observation is important for understanding the physics of neutrino flavor 
conversion in the Sun, since the latest solar neutrino data seems to prefer large mixing 
angles. Here we review the basic arguments and in particular show that the resonance 
does not in general coincide with either the point of maximal violation of adiabaticity in 
the nonadiabatic case or the point of maximal flavor conversion in the adiabatic case. We 
also discuss the modified adiabaticity criterion. 
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1 Introduction 



The MSW mechanism ^, |^ of neutrino flavor oscillations is a simple and very beautiful 
idea. The effect arises because of the coherent forward scattering of the neutrinos on electrons 
and nucleons in matter. This scattering makes the "neutrino index of refraction" of the media 
slightly different from one. The remarkable fact is that this small (and otherwise unobservable) 
change of the index of refraction can have a dramatic effect on the oscillation probability of solar 
neutrinos. In particular, when the flavor off-diagonal terms in the neutrino Hamiltonian are 
small compared to the diagonal terms in vacuum, the flavor conversion can be almost complete. 

It is this possibility of the almost complete conversion that initially attracted a lot of atten- 
tion. The phenomenon was understood as a resonance effect. The flavor conversion was shown 
to occur at the depth in the Sun where the density is such that the matter effects cancel the 
difference between the diagonal terms in the neutrino Hamiltonian. 

The recent data 0, |^, ||, 0], however, together with the latest solar model f^, indicate 
that the preferred value of the mixing angle may be large. It was recently emphasized 
that the resonance picture breaks down in this case. It was shown that the important physical 
processes, level jumping or flavor conversion, do not occur at the resonance. The correct picture 
is to consider adiabatic or nonadiabatic neutrino evolution which reduces to the resonance 
description in the limit 0^1. It is also worth pointing out that this description is perfectly 
continuous across the point of the maximal mixing 6' = 7r/4. 

In this note we review these developments. The aim is to keep the presentation as simple 
as possible, in accordance with the stated goals of the Workshop to make the Proceedings into 
a pedagogical resource. In Section ^ we present an elementary introduction to the resonance 
phenomenon. In Section ^ we pose several questions that seem to be difficult to answer within 
the resonance picture. In Section ^ present the resolutions. 



2 Neutrino oscillations and the MSW effect 

Let us begin by summarizing the standard treatment of neutrino oscillations. This introduction 
is not meant to be exhaustive. In fact, quite often the best presentations of any idea are written 
by people who discovered it or contributed early on, and the case at hand is no exception, see. 



e.g., 0,0, 12, 13, 14]. 



2.1 Simple vacuum oscillations 

Suppose at some point one creates the electron neutrino state, for example as a result of a 
nuclear reaction. Suppose further that there is a term in the Hamiltonian coupling the electron 
neutrino to some other neutrino state. Then, as the neutrino propagates in vacuum, it will be 
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(partially) converted into that state. Still later, the state will again be v^. and so on, i. e., the 
flavor oscillations will take place. 

The simplest case is the case of two flavors. The oscillations occur between the electron 
neutrino and another active neutrino stat^, conventionally labeled v^^. The Schrodinger 
equation for this system is 







= cos6' z/g) — 


sin Q\ 


^2) 


= sin6' z/g) + 


cos 



The presence of the off-diagonal coupling Hep. = H^e ensures that the mass eigenstates \i^i,2) of 
this system are not the same as the flavor eigenstates \i^e,fj)- 

^.)- (2) 

Thus, introducing the off-diagonal term in Eq. (|1|) is equivalent to saying that the mass and 
flavor eigenstates in the lepton sector are not aligned. We know this is the case in the quark 
section, where the two bases are related by the CKM matrix, hence it is not unreasonable to 
expect the misalignment also in the lepton sector. 

The coupling He^ leads to flavor oscillations as the neutrino state propagates in space. If 
at x = the state was purely Ue, after propagating in vacuum to a; = L it will be measured as 
Ufj^ with the probability 

P(z/e ^ /i^) = sin^ 29 sin^ (^^l) . (3) 



Here A is one half of the mass splitting, A = ^JH^^ + {H^^ — ifee)^/4 = Am'^/2E,y. Clearly, 
A sin 2^ = Hefj^, A cos 26' = (i?^^ — Hee)/2. Since all oscillation phenomena depend only on 
(if^^ — ifee)/2 and not on if^^ or Hee separately, we shall often drop the constant term in the 
Hamiltonian 




-A cos 29 A sin 29 
A sin 29 A cos 29 



(4) 



"'^One should keep in mind that the subsequent analysis does not distinguish between i^^ and Vr- Hence the 
notation j/^ does not literally mean the muon neutrino. In fact, the popular oscillation scenario which seems 
to be consistent with the atmospheric and solar neutrino data involves oscillations between and a roughly 



50-50 mixture of and Vt 
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2.2 The resonance and the MSW effect 



2.2.1 Matter effects 

The above picture describes the evolution of the neutrino state in vacuum. One might be 
tempted to also use the same formulas for the evolution in matter. Indeed, it might seem that 
the presence of matter has no effect on the neutrino evolution. After all, the scattering cross 
section for the electron neutrino a ~ GpE^^ is so small for solar neutrinos (-Eiab ~ (0.2 — 15) 
MeV), that the vast majority of them go right through the Sun undeflected. 

Nevertheless, the presence of matter does affect the neutrino evolution in an important way. 
Since the medium is transparent to neutrinos, it is helpful to think in the language of "neutrino 
optics" . The presence of the medium can lead to a index of refraction for the neutrinos. Since 
the index of refraction is a phase phenomenon - it is related to the amplitude of the coherent 
forward scattering - the effect involves the first power of Gp- Wolfenstein pointed out [0 
that this coherent forward scattering leads to the addition of the extra terms to the diagonal 
elements of the Hamiltonian, \/2Gp{Nf, — Nn/2) for the electron neutrinos and \/2Gp{—Nn/2) 
for the other active species. Thus, 

/ A-Acos20 A sin 20 \ 
= const + Asin2^^ Acos2^^-Aj 

where A = GpNe/ \/2. The effect of the neutrino interaction with matter is to modify the 
mixing angle and the mass splitting in matter, 

_ , / A„cos26'm Am sin 29. m \ 
ii^^t - const + A™,cos2^^™, • ^""^ 



where 



Am = - 2AA COS20 + A2, (7) 

A sin 20 



tan 29m = 



A cos 29 - A' 



To find the oscillation probability in a medium of constant density, one can use Eq. (j^) with 
A Am, 9 9m- This results, for example, in the suppression of the oscillations in the solar 
core for A < lO'^ eVVMeV IJ. 



2.2.2 Resonance: the idea 

The true importance of the matter effects for addressing the solar neutrino problem was not 
realized until 1985, when the so-called resonance phenomenon was discovered P, 0]. It was 
known at the time that the Chlorine and Gallium data showed a 50-70% suppression of the 
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R/Rq 

Figure 1: Electron density in the Sun according to the BP2000 solar model. The dashed 
line shows the best-fit exponential ne/NA = 245exp(— lO.54i?/i?0). 



solar neutrino flux. One way to get this suppression was to have a large mixing angle and 
oscillations in vacuum. The resonance idea provided another way. It showed that a large flavor 
conversion could occur even if in vacuum |i^e^t| ^ I-^m/^ ^ -^eel- 

In a nutshell, the idea is to have the matter term cancel the difference of the diagonal 
elements at some depth in the Sun. If, at some depth x, \A[x) — Acos2^| < A sin 26', the 
conversion will take place according to idip^/dt ~ Asin2^(/)e. Notice that the condition 

A(xo) = Acos2^ (9) 

a priori does not require any large fine-tuning of the neutrino parameters, since the electron 
density in the Sun smoothly varies over several orders of magnitude (see Figure [^). 

To get enough conversion, the distance 5x over which the cancellation persists must be long 
enough so that a significant fraction of the oscillation cycle can be completed, 

Asin2^(5x) > 1. (10) 

If one assumes that the profile around the resonance x = xq is approximately linear A{x) ~ 
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^(^^o) + A'{xq)Sx, the width of the resonant region can be written as 



and Eq. ( ]T0| ) becomes 



2Asin2^ , , 

26x ~ . 11 



2A2sin2 20 ^ , , 

> 1. (12) 



A'ixo) 

Whenever the condition in Eq. ([T^ ) is satisfied, one can expect the conversion to be large. 



2.2.3 Resonance: level crossing 

On closer inspection, it turns out that the conversion is not only large, but can be almost 
complete. This can most easily understood as the case of adiabatic level crossing, a phenomenon 
well known in quantum mechanics. If parameters in the Hamiltonian vary slowly enough, the 
state that was initially created as the eigenstate of the instantaneous Hamiltonian (henceforth, 
"the matter mass eigenstate" ) will follow the changing eigenstate. In the case of solar neutrinos 
with a small mixing angle and A < 10~^ eV^/MeV, the fiavor eigenstate i>e in the core almost 
coincides with the heavy matter mass eigenstate 1^2 . The heavy mass eigenstate is, in turn, 
smoothly "connected" to the heavy eigenstate in vacuum which is predominantly composed of 
ly^. The situation is illustrated in Fig. 0. 

To see explicitly how the adiabatic evolution works, let us transform the evolution equation, 
Eq. (|l]), from the fiavor basis (fg, ^fj.) to the basis of the matter mass eigenstates (z/i, z/2), 

= 1/V- 

Remembering that the rotation matrix 



^^/cos^„ -sine \ 
' smfc'm cost/^ ' 



is position dependent, we get 

d 



dx 



dib + dV'^ 

= VHV^-iV^i). 14 
dx dx 



Since VHV^ = diag(-A„, +A„) and from Eq. (0) 

dV^ ^( l\dern 
dx [ -I dx' 
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Figure 2: The level crossing in the Sun (small mixing angle). 



we obtain the desired evolution equation in the basis of the matter mass eigenstates [p^ , 




dOm/dx —iAm 




(15) 



This shows that the coupling between the matter mass eigenstates depends on dOm/dx and, 
hence, vanishes when the density change along the neutrino trajectory is very gradual. More 
precisely, the relevant quantity is the ratio Am/ {dOm/ dx), which at the resonance point equals 
2A? sin^ 26 /A'{xq). Thus the condition for adiabaticity A^/ {dOm/dx) ^ 1 is precisely the same 
as Eq. ([T2D , as should be expected since they both describe the same physics. 



2.2.4 Vacuum oscillation solutions 

If matter effects are important in the Sun, what about the vacuum oscillation solutions to the 
solar neutrino problem? Wouldn't the matter effects also modify the neutrino evolution in that 
case? 

This question can be answered in several steps. First, the vacuum oscillation solutions are 
characterized by the values of Am? in the range 10^^^ — 10^^ eV^. For Am? in this range, the 
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matter effects are obviously important in the core: the term A{x) completely dominates in Eq. 
(|^), leading to the complete suppression of the oscillations there. 

Next, consider what happens in the outer layers on the Sun. If A/Ei, ~ 10~^^ eV^/MeV, 
the off-diagonal term in the Hamiltonian is so small that the neutrino does not have any time 
to oscillate into z/^ while still inside the Sun. In other words, the condition (|12D in this case is 
badly violated and the neutrino evolution is extremely nonadiabatic (flavor conserving). In the 
language of the matter mass eigenstates, this corresponds to the probability of level jumping 
Pc = cos^ 6. One then expects the vacuum oscillations to start up near the surface of the Sun^ 

Finally, it is worth mentioning that this conclusion does not extend to A/E^, ~ 10^^ 
eV^/MeV. A careful analysis shows that the neutrino evolution in this range is not extremely 
nonadiabatic, unlike previously thought, i. e., Pc < cos'^6. For details, see WE. 



3 Questions 

Above we have described several key features of the neutrino evolution in the Sun. They form 
the well established physical picture of the MSW mechanism as the "resonant amplification 
phenomenon" . This picture has become the standard lore over the years. It can be summarized 
as follows: 

• If the evolution is adiabatic, the flavor content of the neutrino state changes as the neutrino 
traverses the Sun. The flavor conversion occurs around the resonance layer, see Eq. (P). 

• To check whether the evolution is adiabatic or not, one should evaluate the adiabaticity 
condition. This condition becomes the most critical at the resonance, Eq. (|12[). If the 
evolution is nonadiabatic, the jumping between the matter mass eigenstates takes place. 
This jumping also happens mostly in the resonance layer. 

This picture is remarkably simple and at the same time very successful in describing the 
neutrino evolution when the mixing angle is small. It, however, encounters serious difficul- 
ties explaining physics in other cases. We shall next pose several questions that expose the 
limitations of this picture. 

3.1 The relevant part of the density profile 

In the preceding discussion we assumed that the density profile was linear. It may not be imme- 
diately obvious why this assumption is justified. After all, the electron density throughout the 
Sun, generally speaking, is not linear. As Figure |l| shows, it falls off more or less exponentially 

^ A detailed analysis of the oscillation phase was recently done in . The value of the oscillation phase in 
the vacuum oscillation region was found to be as if the vacuum oscillations started up at r = O.87i?0. 
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in the radiative zone, while deviating from the exponential near the solar edge and also in the 
core. 

More precisely, we assumed that the density can be well approximated by a linear function 
in some neighborhood (a;o — Sx,xq + 6x) of the resonance point, where 6x = 2A sin 2^/y4'(xo) 
is the previously introduced width of the resonance region. Obviously, this works well when 
A"{xo){6x)y2 < A'{xo){6x), or 

Asm29 <^{A\xo)Y/A"{xo). 

We can see that, when the mixing angle 6 is small, the width of the resonance region is also 
small, so that the linear approximation for the density should be good. The important point is 
that all the "interesting" physics, i. e. the flavor conversion or level jumping, happens within 
this region. The jumping probability in this case can be found using the analytical expression 
for the linear density profile, 

, , / A2sin2 2^\ , , 

(Pc)(iin) = exp I -TT ^,|^_^ 1 ■ (16) 

(Notice that the expression in the exponent is up to a factor just the adiabaticity parameter.) 
The details of the profile outside this region are not so important, so long as the density there 
is sufficiently different from resonant. 

The question is: what happens when the mixing angle is large? According to Eq (|TTp, the 
resonance region in that case is wide and so the linear approximation may no longer apply. 
While the analytical expression for the infinite exponential profile [|r^, [T^, 



g47rroA cos'^ 6 

is not limited to the case of small 6, it is not clear a priori if it can be applied for solar 
neutrinos, since there we are not dealing with an infinite exponential. If the "interesting" 
physics, in the sense defined before, happens on the part of the trajectory which does not lie 



within the exponential part of the profile, we have no justification to use Eq. i\L7\) . 

The question is not just academic, but has a direct application for the present-day solar 
neutrino analysis. Recent solar neutrino fits indicate that the so-called Large Mixing Angle 
(LMA) solution, characterized by Am^ ~ 10~^ — 10^^ eV^ and tan^ ~ 0.1 — 0.9, gives 
a marginally better fit than the small mixing angle (SMA) solution. As the mixing angle 
approaches 7r/4, the resonance defined in Eq. (|]) occurs closer and closer to the solar surface. 
Does this mean we need to know the details of the profile in the convective zone of the Sun 
to reliably compute the neutrino survival probability for the LMA solution? The resonance 
amplification picture says "yes" , while the numerical computations say "no" . 
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3.2 The case of "inverted hierarchy" 

The next issue is understanding the evolution in the case A < ("inverted hierarchy"). This 
means the heavy mass eigenstate in vacuum is made up predominantly of Ve- Since in this case 
no resonance occurs, it received comparatively little attention in the literature until recently. 
At the same time, this case raises several interesting questions that are not easy to answer 
within the traditional resonance framework described above. 

First of all, notice that the mixing angles at the production point in the core and in vacuum 
are different, just like in the A > case. This means that if the evolution is adiabatic, the 
flavor composition of the neutrino state changes as it travels through the Sun. Where does this 
change occur if there is no resonance? 

Next, consider the limit of small A. For A ~ 10"" - 10"^° eVVMeV one expects vacuum 
oscillations which are no different then in the A > case. Namely, the flavor oscillations in 
the solar core are suppressed and in the region close to the solar surface the off-diagonal terms 
in the Hamiltonian are too small to allow a significant fiavor conversion within the Sun. 
Hence, viewed in the mass basis, the evolution in this case is extremely nonadiabatic. The 
jumping between the matter mass eigenstates does occur, even if there is no resonance. How 
can this be? 

4 Matter effects and large mixing angles 

4.1 Connecting A < and A > parts of the parameter space 

The two cases, A < and A > 0, the way they were described up to this point, might 
appear qualitatively very different. Nevertheless, it can be easily demonstrated that they are 
continuously connected. 

To see this, notice that the situations 1)A<0, ^^<7r/4 and 2)A>0, ^^>7r/4 are physically 
the same [|19|, |2^. This can be easily seen from Eqs. (|1]) and (^). The case Q ^ -n jl — Q reduces 
to A ^ —A by redefining the sign of, say, 0^. Hence the physical parameter space can be 
continuously parameterized by keeping the sign of A fixed and varying Q from to 7r/2. This 
means that there should be a unifying physical description of the MSW effect that continuously 
incorporates the resonant {6 < 7r/4) and the nonresonant {6 > vr/4) cases. 

4.2 The points of maximal violation of adiabaticity and maximal 
flavor conversion do not coincide with the resonance 

Such description is not difficult to construct. The two questions we must be able to answer are: 

• In the adiabatic case, what defines the center of the region where the fiavor conversion 
takes place? 
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Figure 3: Neutrino state evolution in the case of the infinite exponential density profile for 
IS.W? jEy = 10~^ eV^/MeV. The plot shows the probability of finding the neutrino in the heavy 
matter mass eigenstate 1^2 as a function of position x, for four different values of the vacuum 
mixing angle. The points of the maximal violation of adiabaticity, as predicted by Eq. (pSf), 
are marked. 



• In the nonadiabatic case, what defines the center of the region where the jumping between 
the mass eigenstates occurs? 

Once again, the description we are seeking should be continuous across the maximal mixing 
[6 = 71/4) and should reduce to the standard resonance picture in the small 6 limit. 

The answers to both questions turn out to be straightforward. First, consider the adiabatic 
case. In this regime the neutrino state is "glued" to the changing mass eigenstate. Hence, one 
can define the point a which the mass basis rotates at the maximal rate with respect to the 
flavor basis to be the point where the flavor composition of the state changes at the fastest 
rate. In other words, we need the point where 6m = dOm/dx is maximal. 

To get some understanding of where this point lies, let us flnd it explicitly for the inflnite 
exponential proflle A{x) oc exp(— x/rg). Using A'{x) = —A{x)/rQ and 

• _ sin^26'^ dA 

~ 2Asin2^t/x' ^ ' 
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^ A sin(2g^^ - 26) 

we find that 9m oc sin2^m sin(2^m — 2^^) = (cos 20 — cos(40m — 20))/2, and so the maximum 
occurs when 

0-- = 7r/4 + 0/2, (20) 

i e., when 0.^ is halfway between its value deep in the core (7r/2) and its value in vacuum {9). 
The corresponding value of the matter term at this point is 

A{9Zn = A. (21) 



Notice that conditions in Eqs. (|20D and (|2lD are completely general, valid for all values of 
the mixing angle. Notice also that they agree with the resonance description 

A{xo) = Acos29, (22) 
9mi^o) = vr/4, (23) 

in the limit ^ 1, as expected. Thus, the resonance amplification interpretation is just a 
particular limit of this general physical picture. 

Now consider the nonadiabatic case. The "level jumping" region will be centered around 
the point of maximal violation of adiabaticity, which corresponds to the minimum of Am/9m, as 



discussed in Sect. |2.2.3| . Since 9^ ^ oo deep inside the Sun and also in vacuum, the minimum 
of Am/9m, unlike the resonance, exists for all values of 9. 

Once again, let us demonstrate this on the example of the exponential profile. A short 
calculation shows that the minimum in this case occurs when 

cot(20„, -29) +2 cot{29m) = 0, (24) 

or 

A = A£2£^^±i8±£°5!2?. (25) 
4 

Thus, the density at the center of the nonadiabatic part of the neutrino trajectory varies between 
^4 = A (0 — > 0) and A = A/2 [9 vr/2)0. The resonance description again only works when 
9 <^ 1. As0— i>7r/4, the resonance moves to infinity, while the true point of maximal violation 
of adiabaticity moves only to the point where A = A/\/2. 

^This is simple to understand. The term l/9m has a minimum aX A = A, and since is falling with x, 
the minimum of the product in the case of large 6 is shifted towards somewhat smaller A. When 6 is small, the 
minimum of l/9m is very sharp and completely determines the minimum of the product. 
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The situation is illustrated in Fig. which shows the probability of finding the neutrino 
in the heavy mass state 1^2 as a function of the distance x. The parameters of the exponential 
were taken from the fit line in Fig. and Arn^ / = 10^^ eV^/MeV. Three large values of 
the mixing angle {6 = 7r/6,7r/4, and vr/3) and one small value {6 = 7r/60) were chosen. The 
dashed lines and dots mark the points where adiabaticity is maximally violated, as predicted 
by Eq. (p5[). One can see that the partial jumping into the light mass eigenstate in all four 
cases indeed occurs around the marked points. 



4.3 Adiabaticity criterion for large mixing angles 

At last, we shall formulate the adiabaticity criterion that is valid for all values of the mixing 



angle. From the preceding discussion it should be obvious that the criterion of Eq. (12), 

\Am/9m\e^=^/i > 1, (26) 

should only work in the hmit ^ ^ 1. At the very least, the quantity Am/Om should be evaluated 
at the true point of the maximal violation of adiabaticity. As an estimate, we can approximate 
it by the point halfway between 9 and 7r/2 (see Eq. (|20D). We then get 



\A,n/e > 1. (27) 

This is better, but still not quite right for large 6. The key is to express the information 
contained in the system of two evolution equations in a single equation. Eq. (|I]) contains 
redundant information: while it contains four real functions, the actual number of degrees 
of freedom is only two, since the normalization of the state is fixed and the overall phase 
is physically irrelevant. It will be easier to formulate the criterion in question once the two 
physical degrees of freedom are isolated. 

With this in mind, let us write down the evolution equation for the ratio s = From 
Eq. (^) we obtain the following first order equation: 

ds 

— = 2tAr^s-em{s^ + l), (28) 

or 

^ = 2^^s-{s' + l). (29) 

dUm Orn 

It is easy to see that the adiabatic limit corresponds to neglecting the terms in parentheses, 
while the extreme nonadiabatic limit is obtained if one neglects the first term on the right. In 
the second case the solution is s = cot{6m)- Thus, the self consistent condition to have the 
extreme nonadiabatic solution is 2\Am/0ms\ <^ (s^ + 1), with s = cot{9m), 6'm = 7r/4 + 9/2. In 
the opposite limit, the evolution is adiabatic. The adiabaticity condition is then 

\AjL\e^^n/4+e/2 > (tan(7r/4 + 6/2) + cot(7r/4 + e/2))/2, (30) 
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Figure 4: The boundary between the adiabatic and nonadiabatic regions in the case of the 
exponential profile, as predicted by Eq. (p6D, Eq. (|27|), and Eq. ([3T|). The dashed curves show 
the contours of constant Pc- For large angles Eq. (^) provides a better description of the 
boundary. 
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or, upon simplification, 

>1. (31) 

Eq. (0) is valid for all values of the mixing angle. 

To illustrate this result, let us again consider the case of the infinite exponential profile. 
Eq. (^) in this case becomes 

8Arosin2 0>l. (32) 

The solid line in Fig. ^ shows the contour of 8Aro sin^ 6 = 1 computed for tq = -R0/IO.54 (the 
fitted value in the BP2000 solar model, see Fig. |ip. For comparison, the dash-dotted curves 
shows the corresponding prediction of Eqs. (p6|) and (^). The dashed curves are the contours 
of constant computed using Eq. ([T7|) . It is clear from the Figure that the description of Eq. 
([3l| ) is correct not only for small 9, but also for 9 ^ 7r/4. 



5 Conclusions 

We have seen that the neutrino evolution in the Sun can be described as the resonance phe- 
nomenon only in the limit 9 <^ 1. For large values of 9, including 9 > 7r/4, the resonance 
condition does not correspond to anything physical and one instead should simply be thinking 
about either adiabatic or nonadiabatic evolution. 

Since the physics is continuous across 9 = n/A, there is no reason to truncate the solar 
neutrino parameter space at maximal mixing. In fact, the latest solar neutrino fits correctly 
present their results in the < 9 < 7i/2 parameter space P^ . 

Notice, that so long as the neutrino survival probability is computed numerically using the 
full solar density profile, the results are unaffected by the considerations presented here. The 
value of the preceding arguments is that they give us a correct physical picture of what is going 
on when neutrinos travel through the Sun and, for example, explain why in the case of the 
LMA solution the results are not sensitive to the details of the profile in the convective zone, 
or why the quasivacuum solutions smoothly extend into the region 9 > tt/A. To paraphrase 
Richard Feynman, we understand the neutrino evolution only if we can predict what should 
happen in various circumstances without actually having to solve the evolution equations each 
time. 

Finally, it is worth mentioning an important article by A. Messiah ||2^ which was brought 



to the attention of the author after Ref. had been submitted to Physical Review. In the 
article, Messiah considers the mixing angles continuously from to it/2 and determines the 
location of the point of maximal violation of adiabaticity for all values of 9. Even though his 
adiabaticity criterion disagrees with the one presented here, it was clearly a very important 
work that, unfortunately, was almost completely forgotten in the subsequent literature. 
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